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ABSTRACT 

A mathematical  stand-outbreak  model  has  been  developed  to  simulate 
population  dynamics  of  the  Douglas-fir  tussock  moth,  the  associated 
defoliationy  and  its  effect  on  host  trees  during  and  immediately  after 
an  outbreak.  The  model  is  available  as  a collection  of  programs  on  the 
USDA  Fort  Collins  Computer  Center  UNIVAC  1100  computer . Information  is 
provided  here  for  use  of  the  model  on  that  computer.  Methods  for 
utilizing  various  options  are  also  discussed. 

I.  INTRODUCTION 

A stand-outbreak  model  has  been^d4vel oped  to  simulate  outbreaks  of 
the  Douglas-fir  tussock  mot'R/fDFTM) , jQrgyia  pseudotsugataT  (McDunnough) , 
a defoliator  of  Douglas-fi^v[£seudotsuga  menziesi jj  var.  glauca  (Beissn.) 
Franco,  and  the  true  fir^I^bTes  spp^,  in  western  North  America.  This 
model  was  developed  by  scientists  funded  by  the  USDA  Expanded  Douglas- 
fir  Tussock  Moth  Research  and  Development  Program.  Details  of  the 
mathematical  form  of  the  model  are  given  by  Colbert,  Overton,  and  White 
(1979a). 

This  model  is  available  at  the  USDA  Fort  Collins  Computer  Center 
(FCCC)  as  a package  of  programs  written  in  ASCII  FORTRAN  for  the  UNIVAC 
1100  system.  The  Forest  Insect  and  Disease  Management  (FI&DM)  Methods 
Application  Group  (MAG)  in  Davis,  California,  has  been  assigned  respon- 
sibility for  maintaining  this  and  other  programs  as  a national  system 
and  to  provide  training  and  documentation  for  potential  users. 


1 This  manual  was  partially  developed  under  the  support  of  the  USDA 
Expanded  Douglas-fir  Tussock  Moth  Research  and  Development  Program, 
FS-PNW  grant  no.  51. 

2 Research  Coordinator,  CANUSA  Spruce  Budworm  Program-West , USDA 
Forest  Service,  PNW  Station,  Portland,  Oregon. 

3 Mathematical  Statistician,  USDA  Forest  Service,  Forest  Insect 
and  Disease  Management,  Methods  Application  Group,  Davis, 

Cal  i form' a. 


This  report  is  designed  as  a computer  user's  manual.  It  is  assumed 
that  those  reading  it  have  some  familiarity  with  computers;  those  inter- 
ested in  the  field  procedures  as  applicable  to  this  model  should  refer 
to  Colbert  and  Campbell  (1979). 

A set  of  procedures  is  provided  in  this  manual  to  illustrate  how  to 
use  the  model.  There  are  three  steps  in  the  use  of  the  stand-outbreak 
model : 

-Decide  which  program  options  are  desired. 

-Prepare  input  files  that  supply  information  to  the  program(s). 

-Develop  an  appropriate  set  of  instructions  to  the  computer  (called 
a "runstream") . 

A discussion  or  the  development  and  present  structure  of  the  stand- 
outbreak  model  follows  to  assist  in  understanding  the  options  available 
and  their  use. 

A.  Model  Structure 


The  model  is  organized  into  two  levels  (Fig.  1).  At  the  top  level 
is  the  stand-outbreak  module  (Sq),  which  has  the  least  detail.  At 
this  level,  spatial  resolution^  is  the  entire  stand,  and  temporal 
resolution  is  annual  and  spans  up  to  10  years,  including  four  years  of 
insect/host  interaction.  In  the  model,  an  outbreak  episode  consists  of 
five  phases.  The  first  three  are  the  three  years  of  the  outbreak 
sequence:  release,  peak,  and  decline  (Wickman  et  al . 1973).  Because 
residual  population  densities  can  still  be  high  at  the  end  of  the  third 
phase,  a post-decline  fourth  phase  has  been  added  to  the  outbreak  model. 
To  assess  the  effects  of  defoliation,  the  stand  is  followed  for  another 
six  years  in  phase  5. 


Host  trees  are  grouped  and  sorted  into  classes.  A tree  class  is 
defined  by  tree  species;  nominal 5 foliage  conditions;  actual  foliage 
conditions  at  the  beginning  of  the  simulation;  and  insect  density  at  the 
beginning  of  the  simulation.  Each  tree  class  becomes  a module  at  the 
second  level  (S^)  in  the  hierarchy.  These  modules  operate  indepen- 
dently of  each  other  at  annual  resolution,  except  that  redi stribution  of 
insects  among  trees  within  the  stand  is  simulated  by  annually  reducing 
the  variance  in  insect  density  among  classes.  The  user  decides  how  many 
tree  classes  to  use  and  how  to  partition  the  stand. 


4 Resolution  refers  to  the  degree  of  detail  that  can  be  observed  at  a 
particular  level  in  the  model  structure. 

5 Throughout  the  manual,  "nominal"  refers  to  conditions  before  any 
defol i at  ion. 
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Figure  1.  Stand-outbreak  model  structure  and  resolution. 


Insect-foliage  interactions,  insect  population  dynamics,  and  foliage 
dynamics  are  followed  through  the  model  branch,  which  operates  at  the 
most  detailed  resolution.  This  is  the  operational  portion  of  the  tree 
class  module.  The  model  branch  is  defined  as  1000  sq.  in.  of  nominal 
midbole  branch  area  along  with  its  insect  population.  This  module 
follows  the  insect  and  foliage  through  each  year  by  successive  projec- 
tions of  new  and  old  foliage  biomass,  and  number  of  insects,  and  larval 
biomass.  The  temporal  resolution  at  this  level  is  the  "occasion".  Each 
year  is  divided  into  10  occasions.  The  first  six  of  these  are  10-day 
intervals  of  instar  development,  feeding,  and  mortality.  Insects  and 
foliage  are  followed  through  the  rest  of  the  year  in  the  remaining  four 
occasions. 

B.  Scope 


The  model  was  developed  using  data  from  an  outbreak  which  occurred 
in  the  Blue  Mountains  of  northeastern  Oregon  during  1971-74.  The  model 
reflects  both  insect  life  table  information  and  observed  effects  of 
insects  on  host  trees.  Host  species  were  Douglas-fir  and  grand  fir,  A. 
grandi s (Dougl.)  Lindl..  At  present  only  limited  information  is  avail- 
able on  other  host  species. 

The  model  can  simulate  any  size  stand  in  as  much  detail  as  is 
available  in  the  data  base.  No  restriction  exists  on  the  number  of  tree 
classes  that  can  be  used.  Default  parameters  cause  the  model  to  simu- 
late the  1971-74  Blue  Mountains  outbreak,  but  minor  parameter  changes 
can  cause  it  to  behave  like  outbreaks  that  have  occurred  elsewhere 
providing  that  data  are  available  (Colbert  et  al . 1979b).  As  tussock 
moth  life  table  information  becomes  available  from  other  areas  and  the 
model  receives  more  use,  its  scope  should  increase.  Each  future  out- 
break will  provide  information  that  will  allow  the  model  to  reproduce 
local  eccentricities  (i.e.,  deviations  from  the  general  pattern  estab- 
lished for  the  Blue  Mountains)  and  further  validate  the  outbreak 
epi sode. 

With  the  advent  of  early  detection  of  outbreaks  provided  by 
pheromone  traps,  complete  outbreak  episodes  can  be  followed.  With  this 
information  and  the  details  of  effects  of  defoliation  on  other  hosts 
such  as  white  fir,  A.  concol or  (Gord.  and  Glend.)  Lindl.,  and  subalpine 
fir,  A.  1 asiocarpa  THook. ) Nutt.,  the  model  could  be  calibrated  to  simu- 
late tussock  moth  outbreaks  accurately  wherever  they  occur. 

C . Appl i cations 

The  stand-outbreak  model  was  originally  developed  to  help  organize 
ecological  research.  It  has  evolved  into  what  is  hoped  will  be  a useful 
tool  for  biological  evaluations  and  forecasts  of  population  trends  and 
their  effects.  Alone,  it  projects  both  a typical  outbreak  and  the 
effects  of  defoliation  in  terms  of  subsequent  growth  loss,  top-kill,  and 
mortality.  By  exercising  some  of  the  options  (Section  V),  it  can  also 
project  the  effects  of  potential  control  al ternatives.  The  effects  of 
variations  in  both  natural  mortality  and  other  environmental  factors  can 
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also  be  explored.  For  example,  changes  in  "normal"  virus  loads  can  be 
accommodated  by  altering  larval  disease,  and  pupal  and  overwinter 
mortality  rates. 

Long-range  effects  of  changes  in  stand  composition  and  stocking  on 
tussock  moth  outbreaks  can  be  explored  when  the  stand-outbreak  model  is 
used  in  conjunction  with  the  stand-prognosis  model  (Campbell  and 
McFadden  1977).  This  system  is  operational  at  FCCC  and  is  maintained  by 
the  USDA  Forest  Service  Timber  Management  Staff  at  Fort  Collins, 

Col orado. 


II.  MODEL  INPUT-OUTPUT  STRUCTURE 

Model  input  is  information  the  computer  must  have  to  complete  a 
simulation.  Four  data  sets  or  files  are  needed.  Two  of  these  deal  with 
defoliation  effects,  one  for  each  host  species.  These  two  files  contain 
information  utilized  by  the  model  for  post-outbreak  evaluations.  The 
remaining  files  are  structured  such  that  initial  insect  densities  and 
foliage  conditions  on  the  model  branch  are  contained  in  one  file  while 
the  other  file  contains  information  on  insect  mortalities  and  other  mod- 
el coefficients.  The  information  contained  in  these  two  files  provides 
the  basis  for  describing  the  interaction  between  the  insect  and  the  host 
throughout  the  four  phases  of  an  outbreak. 

All  information  produced  through  a simulation  is  termed  model 
output.  During  a simulation  three  tables  are  produced:  a summary  of 
annual  insect  and  foliage  changes  on  the  model  branch  (Table  1);  defoli- 
ation of  the  model  branch  and  tree  and  subsequent  effects  (Table  2);  and 
a summary  of  the  parameters  (quantities  which  determine  the  character- 
istics of  an  outbreak)  that  produced  the  simulation  (Table  3).  Most 
details  of  a simulation  are  also  available  as  an  option  (Table  4). 

A.  Information  Needed 


Four  data  files  are  required  to  produce  a simulation.  Three  of 
these  are  parameter  files.  The  fourth  contains  tree-class  parameters 
and  initial  conditions. 

The  first  parameter  file  is  *PARAMETERS. 6.  Table  3 provides  an 
example  of  most  of  the  information  required  for  this  file.  All  control 
mortality  rates  are  entered  here^. 

The  other  two  parameter  files,  DFTM*DFTMGF.  and  DFTM*DFTMDF. , 
contain  values  for  determining  effects  of  defoliation  on  the  two  host 
species  for  which  data  are  available.  All  parameters  in  the  file 
DFTM*DFTMGF.  are  for  determining  the  effects  of  defoliation  on  grand 


6 File  names  at  FCCC  always  end  with  a period. 

7 All  parameters  in  ^PARAMETERS. , their  default  values,  and  a brief 
description  of  each  parameter  are  listed  in  Appendix  A. 
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fir;  DFTM*DFTMDF.  is  for  Douglas-fir  (Figure  9,  Appendix  D).  While 
these  two  files  are  available  to  the  user,  it  is  anticipated  that  most 
users  will  not  need  to  modify  data  contained  in  them. 

Until  more  complete  data  for  the  effects  of  tussock  moth  defoliation 
on  white  and  subalpine  fir  become  available,  you  may  substitute  the 
information  available  for  grand  fir.  Existing  research  on  white  fir 
indicates  close  agreement  between  white  and  grand  fir  regarding  the 
effects  of  tussock  moth  defoliation. 

Limited  data  from  subalpine  fir  suggest  that  the  moth  affects  this 
species  less  than  grand  fir.  The  stress  mortality  rate  for  tussock  moth 
feeding  on  subalpine  fir  is  higher  than  for  any  of  the  other  host  spe- 
cies (Beckwith  1976,  1978).  The  increased  insect  mortality  causes  less 
defoliation  from  the  same  initial  population.  Hence  less  defoliation 
results  in  less  tree  mortality,  top-kill,  and  growth  loss.  To  date, 
this  information  has  not  been  sufficiently  quantified  to  justify 
separate  parameters  for  this  species,  however. 

The  final  set  of  information  needed  to  produce  a simulation  contains 
the  tree-class  parameters  and  initial  conditions.  These  are  given  by 
specifying  nine  items  for  each  tree  class  (Fig.  2).  This  information  is 
stored  in  a file  named  *IC.. 

A subset  of  the  tree  classes  constructed  as  file  *IC.  can  be 
selected  to  be  used  in  a simulation.  Variation  in  foliage  and  insect 
population  data  can  be  used  to  produce  a set  of  tree  classes  that  will 
describe  the  range  of  potential  defoliation. 

B.  Information  Produced 


The  model  simulation  produces  three  tables,  the  first  two  structured 
by  tree  class.  The  pages  of  Tables  1 and  2 are  numbered  1.1,  1.2,...; 
2.1.  2.2,...  etc.,  to  assist  the  user  in  cross-referenci ng  between 
insects,  defoliation,  and  the  results  of  defoliation.  Table  3 is  an 
example  of  the  summary  of  parameters  from  *PARAMETERS. . 

Table  1 is  an  annual  resolution  summary  of  the  four  outbreak  phases 
(years).  Each  phase  is  summarized  by  displaying  the  number  of  insects 
and  foliage  complement  at  the  beginning  of  the  year,  and  the  defoliation 
and  number  of  viable  eggs  produced  at  the  end  of  the  year. 

Table  2 is  a defoliation  summary  and  effects  table  that  gives,  for 
each  tree  class,  the  giodel  branch  defoliation  for  each  year  of  the 
simulation  and  the  translation  of  maximum  branch  defoliation  into  per- 
cent of  crown  totally  defoliated.  Crown  defoliation  is  then  used  to 
produce  expectations  of  direct  and  secondary  mortality,  top-kill,  and 
growth  loss  to  trees  of  each  class. 

Direct  mortality  is  that  which  is  attributable  to  defoliation  by  the 
tussock  moth.  Secondary  mortality  is  that  attributable  to  bark  beetles 
and  all  other  causes.  Secondary  mortality  is  computed  as  a function  of 
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Figure  2.  Intial  conditions  and  parameters  used  to  describe  a tree 
class  (information  stored  in  the  file  *IC.). 

Item  No.  Description 

1.  Species  code:  1 for  Douglas-fir,  2 for  grand  fir. 

2.  Number  of  trees  or  trees/acre  that  this  tree  class  represents.  Used 
when  redistribution  of  insects  between  tree  classes  is  considered 
(Section  V-D). 

3.  Weight  factor.  Whole-tree  crown-foliage  weight  (kg)  for  one  tree 
representative  of  the  class.  Used  when  redistribution  of  insects 
betv/een  tree  classes  is  considered  (Section  V-D). 


4.  Nominal  percent  new  foliage.  The  percent  of  total  foliage  biomass 
(needle  dry  weight,  g/ 1 000  sq.  in.)  that  is  all  current  year's 
growth.  Nominal  percentage  is  defined  as  current  year's  growth  from 
midcrown  samples  of  branches  that  have  not  been  fed  upon  (Section 
V-D). 

5.  Nominal  total  foliage  biomass  (needle  dry  weight,  g/ 1000  sq.  in.)  of 
midcrown  sample  branches  (Section  V-D). 

6.  Actual  new  foliage  (current  year's  growth)  biomass  at  the  initiation 
of  simulation  (needle  dry  weight,  g/ 1000  sq.  in.).  If  some  previous 
defoliation  is  assumed  or  has  been  measured,  this  will  be  less  than 
the  nominal  amount  from  items  4 and  5 above  (Section  V-E). 


7.  Actual  old  foliage  biomass  at  the  initiation  of  simulation  (needle 
dry  weight,  g/ 1 000  sq.  in.)  (Section  V-E). 

8.  Initial  number  of  established  first  instars  per  model  branch  (1000 
sq.  in.)  (Section  V-F). 


9. 


Tree-class  number.  A reference  or  index 
in  classification.  It  is  used  in  output 
referencing  (Section  V-A). 


number  to  assist  the  user 
tables  for  cross- 
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the  expected  proportion  of  trees  in  each  of  the  five  top-kill  classes 
subsequent  to  direct  mortality.  Once  secondary  mortality  has  been 
removed,  growth  loss  is  computed  for  the  remaining  trees.  Currently  in 
the  model  height-growth  loss  is  proportional  to  radial -growth  loss. 

Table  3 contains  a summary  of  important  parameters  used  to  produce 
Tables  1 and  2 discussed  above.  They  are  the  insect  natural  and  control 
mortality  rates  for  each  life  stage  and  year  of  the  simulation;  nominal 
(maximum)  fecundity;  mean  daily  individual  larval  growth  rate;  the 
destruction/consumption  ratio  (used  to  determine  the  amount  of  foliage 
actually  destroyed  as  the  insects  consume  foliage);  and  the  redistri- 
bution coefficient  and  mean  fall  egg  densities  over  the  stand  for  each 
phase.  The  full  set  of  parameters  used  in  a simulation  are  given  in 
Section  III. 

C.  Model  Details 


During  a simulation,  four  variables  (new  and  old  foliage  biomass, 
number  of  insects,  and  mean  larval  biomass)  are  followed  through  10 
occasions  or  time  intervals.  Values  of  these  variables  at  points 
between  these  occasions  can  be  obtained  by  an  optional  sorting  program. 
This  program  produces  the  additional  table  (Table  4)  of  requested 
details  of  a simulation.  One  variable,  mean  larval  biomass,  was  not 
presented  in  the  earlier  tables.  This  variable  follows  a preset  tra- 
jectory (course)  according  to  the  growth  rates  given  in  Table  3,  as  long 
as  foliage  is  available  to  maintain  growth. 

Figure  3 describes  the  occasion  structure  and  times  at  which 
detailed  information  is  available  for  Table  4.  Occasion  as  used  here 
refers  to  an  interval  of  time  during  which  specific  prescribed  activ- 
ities take  place.  By  stipulation  of  the  variable  of  interest  and  the 
time  during  the  simulated  year  of  interest,  the  Table  4 program  will 
sort  through  the  simulation  details  and  assemble  the  requested  informa- 
tion for  each  tree  class. 

III.  INPUT  FILE  FORMATS 


A.  Parameter  File 


The  parameter  file  is  read  from  the  file  named  *PARAMETERS.  in  the 
following  formats:  first  record  (112,  114,  2013);  remaining  24  records 
(6F12.7).  Figure  4 gives  default  values  for  this  file  in  these  formats. 
The  information  contained  in  this  figure  was  utilized  to  generate  Tables 
1,  2,  and  3 discussed  earlier.  Appendix  A gives  location,  value,  and 
descriptions  for  each  parameter. 

B.  Host  Effects  Files 


The  two  host  effects  files  are  DFTM*DFTMDF.  and  DFTM*DFTMGF.  for 
Douglas-fir  and  grand  fir  respectively.  The  two  standard  format  default 
files  are  given  in  Appendix  D. 
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Figure  3.  The  occasion  structure  of  one  year  of  a simulation  and  the 
times  at  which  detailed  information  is  available. 

Access  times*  Occasion  descriptions 


(initiation,  yr  j)0 


bud  burst,  new  shoot  elongation,  egg 
hatch,  establishment  of  first  instars 


1 

first  instars  feeding,  growth,  and 
mortality 

2 

second  instars  feeding,  growth,  and 
mortal ity 


third  instars  feeding,  growth,  and 
mortal ity 


4- 


& 


& 


7- 


' & 


9 


(initiation,  yr  j+1 ) 0 


fourth  instars  feeding,  growth,  and 
mortality 


fifth  instars  feeding,  growth,  and 
mortality  (after  all  feeding  and 
growth,  males  pupate  and  are  no  longer 
expl icitly  modeled) 


sixth  instars  (females  only)  feeding, 
growth,  and  mortality 


female  pupation,  pupal  mortality, 
adult  emergence,  and  adult  mortality 


egg  laying,  aging  of  foliage,  and  new 
foliage  potential  set 


overwinter  mortality  of  eggs,  shedding 
of  old  foliage,  and  redistribution  of 
insects  between  tree  classes 


The  value  of  the  access  time  given  here  is  one  less  than  the 
occasion  code  for  accessing  the  state  variables  (see  Fig.  8)  (see 
Section  V-I). 
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Figure  4.  Default  values  for  *PARAMETERS.  given  in  record  image 
format. 
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The  other  host  species  can  be  accommodated  by  substituting  any  white 
or  subalpine  fir  mortality  information  into  a file  properly  formatted 
and  used  in  place  of  either  DFTM*DFTMDF.  or  DFTM*DFTMGF.  or  both.  If 
this  is  to  be  done,  then  the  appropriate  species  code  in  the  initial 
condition  file  will  make  the  appropriate  substitutions.  This  will  also 
cause  the  species  code  selection  of  other  species-specific  parameters 
during  simulation. 

C.  Initial  Conditions  File 


The  initial  condition  file,  named  *IC.,  is  read  in  the  format  (112, 
7F9.5,  113).  Figure  2 gives  descriptions  of  entries  of  any  record  of 
this  file.  This  file  is  read  from  logical  unit  1,  4,  or  7 depending 
upon  whether  the  simulation  starts  with  phase  I,  II,  or  III.  An  example 
of  an  initial  condition  file  for  a simulation  starting  with  phase  I and 
covering  the  range  from  1 to  15  initial  insects  per  1000  sq.  in.  is 
given  in  Figure  5.  Each  tree  class  in  this  example  is  Douglas-fir. 
Fields  2 and  3 have  zero  entries  for  all  tree  class  records  because 
insect  redistribution  is  not  considered.  The  model  branch  for  each  tree 
class  is  assumed  to  have  200  g.  of  foliage,  of  which  25%  is  current 
year's  growth  and  with  no  previous  defoliation.  This  file  was  utilized 
in  conjunction  with  the  parameters  file  in  Figure  4 to  produce  Tables  1, 
2,  and  3. 

A sorting  routine  has  been  built  into  the  model  program  such  that 
whenever  any  one  of  the  two  nominal  foliage  input  variables  in  *IC. 
changes  values,  output  control  is  automatically  directed  to  a new  page. 
The  user  should  therefore  group  the  input  records  accordingly  to  elimi- 
nate extraneous  paginations. 

D.  Table  4 

The  data  needed  to  generate  Table  4 are  contained  in  three  files: 
*LUN10. , *LUN25.,  and  *SPECTAB4..  Specifications  of  information  desired 
(Section  V-D)  are  put  in  *SPECTAB4..  Figure  6 is  an  example  of  the 
*SPECTAB4.  file.  This  file  was  used  to  generate  Table  47  The  first 
record  or  card  is  in  format  (113,  112),  and  each  succeeding  record  is  in 
format  (212)  *LUN10.  and  *LUN25.  are  computer-generated  during  a simu- 

lation; *SPECTAB4.  is  generated  by  the  model  user. 

IV.  RUNSTREAM  GENERATION 

Figure  7 provides  a description  of  the  general  process  of  producing 
input  files,  executing  a simulation,  and  obtaining  output. 

The  series  of  computer  commands,  known  as  "control  statements",  that 
results  in  execution  of  a simulation  is  called  the  "runsteam".  At  FCCC, 
the  symbol  0 is  always  the  first  (left  hand)  character  in  a control 
statement.  This  symbol  must  appear  in  column  1.  This  section  outlines 
procedures  for  generating  a runstream  and  the  available  runstream 
options.  We  assume  that  you  will  be  communicating  with  FCCC  through  an 
interactive  terminal  and  that  FCCC  has  assigned  you  an  identifier  (site 
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Figure  5.  The  initial  condition  file  *IC.  for  a simulation  starting 
with  phase  I initial  conditions. 
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Figure  6.  The  input  specification  file  *SPECTAB4.  that  determines 
what  is  included  in  Table  4. 
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Figure  7.  General  process  of  producing  input  files,  executing  a 
simulation,  and  obtaining  output. 
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i.d.)  that  allows  you  to  initiate  communications.  As  used  here,  pointed 
brackets  and  the  description  within  them  should  be  replaced  with  appro- 
priate information,  e.g.,  <account>  should  be  replaced  with  your  FCCC 
account  number. 

A.  Run  Initiation  and  Identification 


1.  Initiate  the  run  with  the  following  control  statement: 

@RUN  <runid> ,<account> ,< project id> 

2.  Identify  this  run  and  its  associated  files  with  a "qualifier".  This 
qualifier  will  allow  you  to  identify  the  input  and  output  files  for  each 
run  you  make.  To  use  a qualifier,  you  have  to  know  the  form  of  a file 
name.  File  names  at  FCCC  are  structured  as  follows: 

<qual i fier>*<fil ename>.<elementname> 

One  type  of  file  does  not  use  element  names,  and  is  referred  to  only  as 
<qual if ier>*<f i lenarne>. . Output  produced  by  the  program  is  placed  in 
these  files,  which  are  called  systems  data  formatted  files.  Whenever  a 
file  is  named  (for  example,  by  an  @ASG  command  with  the  appropriate 
options)  without  a qualifier,  an  explicit  or  implied  qualifier  is 
attached  to  the  front  of  the  file  name  equal  to  the  <projectid>.  For 
example,  if  a run  statement  is: 


@RUN 

• 

AB9XXX,  1234567890AB,  NAME 

• 

@ASG , UP 

FILE2. 

(1) 

@ASG  ,T 

FILE17. 

(2) 

(1)  creates  a file  named  NAME*FILE2..  This  file  will  be  catalogued 
(placed  in  a master  file  directory)  so  it  can  be  referenced  after  the 
run  terminates.  This  file  can  be  referred  to  simply  as  FILE2.  in  any 
run  with  <projectid>  = NAME. 

(2)  creates  a file  named  NAME*FILE17. . This  file  is  not  catalogued 
because  of  the  "T"  (temporary)  option  in  the  @ASG,T.  This  temporary 
file  disappears  after  the  run  terminates.  The  file  can  be  referred  to 
simply  as  FILE17..  By  use  of  the  @QUAL  statement,  the  implicit  qual- 
ifier can  be  changed  for  files  referred  to  as  *<f il ename>. . For 
example: 


@RUN  <runid>,<account>,YOU 
0QUAL  ME 
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@ASG,T 

RECORDS. 

(1) 

@ASG  ,T 

*REC0RDS. 

(2) 

@ASG ,UP 

*STUFF. 

(3) 

(1)  creates  a temporary  file  named  Y0U*REC0RDS.  that  can  be  referenced 
later  in  the  same  run  as  simply  RECORDS,  (or  as  Y0U*REC0RDS. ) . 

(2)  creates  a temporary  file  named  ME*RECORDS.  that  can  be  referenced 
later  in  the  run  as  *REC0RDS.  (or  as  ME*RECORDS. ) . The  two  files 
RECORDS,  and  *REC0RDS.  are  completely  separate. 

(3)  creates  a catalogued  (permanent)  file  named  ME*STUFF.  that  can  be 
referenced  later  in  the  run  or  in  another  run  as  ME*STUFF.  (or  simply 
*STUFF.  if  the  statement  OQUAL  ME  is  already  in  the  runstream).  Files 
always  can  be  referenced  by  the  whole  name. 

Currently  available  packaged  runstreams  for  use  with  the 
stand-outbreak  model  (DFTM*RUNSTREAM. PRE1 , DFTM*RUNSTREAM.PRE2, 
DFTM*RUNSTREAM.P0ST1,  etc.)  assume  that  you  have  already  entered  a 
@QUAL<i.d.>  to  identify  the  output  files  for  this  run.  <i.d.>  should  be 
some  word  to  identify  the  run.  This  word  can  be  up  to  12  characters 
long,  using  letters,  numbers,  and  the  symbols  hyphen  (-)  and  dollar  sign 
($).  For  example: 

OQUAL  TESTRUN36-D5 

B.  Parameter  and  Initial  Condition  Files 


Input  files  which  are  already  in  existence  can  be  utilized  to  create 
new  data  files  for  another  simulation.  For  example,  assume  you  have 
previously  run  a simulation  identified  by  the  qualifier  "RUN56"  and  thus 
have  catalogued  files  RUN56*PARAMETERS.  and  RUN56* I C . . You  now  wish  to 
do  a simulation  identified  by  the  qualifier  "RUN57",  and  the  only 
changes  between  the  input  to  these  two  runs  are  in  the  *PARAMETERS. 
file.  Then: 

ORUN  <runid>,<account>,<projectid> 

OQUAL  RUN57 

@ASG ,UP  *PARAMETERS. 

@ASG,UP  *IC. 

@C0PY  RUN56*PARAMETERS. ,*PARAMETERS. 

@C0PY  RUN56*IC. , *IC. 

@ED ,U  *PARAMETERS. 
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[make  changes  to  the *  *PARAMETERS. 
file  using  the  editor  processor] 


EXIT 


You  now  have  a file  RUN 57*1 C . identical  to  RUN 56*1 C . , and  a file 
RUN57*PARAMETERS.  identical  to  RUN56*PARAMETERS. , except  for  the  changes 
made  above. 

C.  Stand-Outbreak  Runstream 


1.  Select  the  appropriate  packaged  runstream  to  ready  the  computer  for 
execution  of  the  simulation.  Select  DFTM*RUNSTREAM. PRE1  for  phase  1 
initiation,  DFTM*RUNSTREAM.PRE2  for  phase  2 initiation, 

DFTM*RUNSTREAM. PRE3  for  phase  3 initiation.  Example: 

@ADD  DFTM*RUNSTREAM. PRE2 

This  command  should  result  in  a string  of  "READY"  messages  appearing  at 
the  terminal.  When  the  prompt  ">"  appears,  the  computer  is  ready  to 
execute  the  simulation  with  phase  2 initiation. 

2.  Execute  the  simulation: 

@XQT  DFTM*OUTBREAK. 

After  execution,  the  printed  output  from  the  simulation  is  stored  in  the 
files  *15.,  *16.,  and  *17.  and  the  Table  4 input  files,  *LUN10.  and 

*LUN25.  have  been  created. 

3.  Select  an  output  option.  Example  1: 

@MSG,N  <50  character  address  line) 

@MSG,N  <50  character  address  line> 


@ADD  DFTM*RUNSTREAM. P0ST1 

The  above  series  of  commands  produces  a printed  listing  of  the  output 
files  *15.,  *16.,  *17.  with  an  address  as  specified  in  the  @MSG,N  com- 
mands. The  files  *15.,  *16.,  *17  are  also  deleted. 
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Example  2: 

PADD  DFTM*RUNSTREAM.P0ST2 
PSYM  *15.,,  <deviceid> 

PSYM  *16.,,  <deviceid> 

PSYM  *17.,,  <deviceid> 

The  command  PADD  DFTM*RUNSTREAM.P0ST2  produces  no  listing  of  the  output 
files,  and  it  is  left  to  you  to  send  the  output  contained  in  *15.,  *16., 
and  *17.  to  an  appropriate  device  (presumably  your  own  line  printer)  via 
the  PSYM  commands.  Files  are  deleted  after  the  PSYM  unless  the  "U" 
option  is  specified.  For  example: 

PSYMjU  *15.,,  <deviceid> 

Multiple  copies  of  the  output  can  be  produced  by  inserting  the 
number  of  copies  between  the  two  commas  in  the  PSYM  command.  For 
example: 

PSYM  *<myf il e>. ,5,<mydevice> 

produces  five  copies  of  *<myfile>.  at  <mydevice> 

4.  Because  *LUN10.  and  *LUN25.  are  catalogued  files,  they  will  be 
available  as  input  to  the  Table  4 program  for  six  days  after  they  are 
created.  After  six  days,  the  files  are  purged  from  the  system  unless 
preserved  via  the  PSAVE  command.  For  example: 

PSAVE  RUN7*LUN10. , 791231 

description  of  contents  of  file  to  be  entered  by  user> 

This  command  will  preserve  the  file  RUN7*LUN10.  until  12-31-79.  The 
file  will  also  carry  with  it  the  description  specified.  This  feature  is 
useful  when  trying  to  identify  files  created  some  time  in  the  past. 

D.  Specifications  for  Table  4 Output 

It  is  assumed  here  that  you  have  completed  a simulation  and,  during 
the  run,  the  files  <qual >*LUN10.  and  <qual >*LUN25.  were  created.  The 
status  of  these  files  can  be  checked,  if  necesary,  by  use  of  the  PSTAT 
command.  For  example,  PSTAT  *DFTM*RUNSTREAM.  results  in  the  following: 

STATUS  2R1  R72-16  >*8/22/79  12:15:23 

QUAL:  DFTM  ACCOUNT:  1105205205  CREATED  : 06/26/79 

FILE:  RUNSTREAM  STATUS  : SAVED  PURGE  : 12/30/80 

CYC  : 1 LAST  REF:  06/26/79 

DESC : DFTM  OUTBREAK  MODEL  RUNSTREAM  FILES 
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If  the  qualifier  on  these  two  files  is  different  from  the  present 
one,  you  should  make  a qualifier  designation  by  inserting  an 
0QUAL  <qual>  into  your  current  runstream.  Next,  create  the  Table  4 
specification  file  described  in  Sections  III -D  and  V-G.  This  file  must 
be  given  the  name  <qual >*SPECTAB4. . This  file  can  be  created  by  use  of 
the  editor  (@ED,I)  or  the  data  command: 

@DATA,I  <filename>. 


DATA  GO  HERE 


PEND 

E.  Table  4 Runstream 


You  are  assumed  to  have  three  files  with  the  same  qualifier: 

<qual >*LUN10. 

<qual >*LUN25. 

<qual >*SPECTAB4. 

The  first  two  were  created  by  a run  of  the  model,  and  the  last  has  been 
created  by  you.  A further  assumption  is  that  you  are  communicating  with 
the  system  via  an  interactive  (demand)  terminal  and  have  set  the 
qualifier  to  match  the  one  on  these  files. 

Two  runstreams  can  be  used  for  the  Table  4 program,  depending  on 
output  requi rements . The  first,  using  DFTM*RUNSTREAM.TABLE4A,  is  for 
use  when  the  output  is  to  be  routed  to  the  FCCC  line  printer  and  mailed 
to  you.  Then,  only  two  types  of  statements  are  necessary: 

@MSG,N  <50  character  (maximum)  address  for...> 

0MSG,N  <... mailing  from  FCCC  to  user> 

0ADD  DFTM*RUNSTREAM.TABLE4A 

This  runstream  will  first  make  all  the  necessary  file  assignments  and 
linkages  and  will  then  execute  the  Table  4 program.  After  execution  of 
the  Table  4 program,  its  output  will  be  disposed  (via  @SYM  commands)  to 
the  FCCC  line  printer  with  an  address  given  in  the  @MSG,N  commands. 

Note  that  this  will  delete  the  output  file  as  it  prints  it.  If  the  user 
desires  to  preserve  the  Table  4 output  file,  he  should  use  the  procedure 
given  below. 
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The  second  procedure,  using  DFTM*RUNSTREAM.TABLE4B,  is  for  use  when 
the  output  is  to  be  routed  directly  to  your  device.  The  runstream  will 
then  be: 

@ADD  DFTM*RUNSTREAM.TABLE4B 

@SYM  *TABLE4. ,<copies> ,<your  device> 

This  runstream  will  make  all  the  necessary  file  assignments  and  linkages 
and  execute  the  sorting  program,  but  disposition  of  the  program's  output 
(contained  in  *TABLE4.)  is  left  to  you.  Here,  too,  if  you  want  to 
retain  the  output  file  generated  (*TABLE4.),  you  must  indicate  the  "U" 
option  on  your  @SYM  commands  and  use  0SAVE  if  you  wish  to  retain  the 
file  for  more  than  six  days. 

V.  OPTIONS  AVAILABLE 

Many  options  are  available.  Most  of  the  relations  that  may  be 
varied  are  assessed  through  changing  the  model  parameters.  Most  of 
these  values  are  read  in  at  the  time  of  execution  and  accessed  through 
the  four  files  described  previously.  Parameter  options  and  access  to 
simulation  details  are  described  and  discussed  in  this  section. 

A.  Years  to  be  Simulated  and  Tree-Class  Selection 


Two  options  are  exercised  through  data  supplied  on  the  first  record 
of  *PARAMETERS.  (Appendix  A):  the  phase  in  which  the  simulated  outbreak 
begins  and  the  number  of  tree  classes  to  be  included  in  the  simulation. 

1.  You  may  start  the  simulation  in  phase  I,  II,  or  III  by  choosing 
the  value  1,  2,  or  3 for  the  first  entry  on  the  first  record  of  the 
parameter  file,  *PARAMETERS.  (Appendix  A;  i =1 , j=l).  If  the  simulation 
is  to  start  in  phase  I,  initial  foliage  complement  is  assumed  to  be  the 
same  as  the  nominal.  Initial  insect  populations  should  be  below  20 
established  first  instars  per  1000  sq.  in.  from  midbole  samples. 

Starting  in  phase  II,  light  defoliation  may  be  visible,  but  no 
heavily  defoliated  patches  should  have  appeared;  model  projections  from 
phase  I (Table  1)  show  only  minimal  deviations  from  nominal  foliage 
conditions.  If  a simulation  is  to  start  in  phase  III,  however,  consid- 
erable previous  defoliation  could  have  occurred.  "Hot  spots"  may  have 
greater  than  90  percent  defoliation  from  the  previous  year.  This  is  the 
year  of  heavy  defoliation;  collapse  is  beginning,  survivorship  is  low, 
and  the  population  trend  from  the  previous  year  is  less  than  1.0. 

2.  The  total  number  of  tree  classes  desired  should  be  put  in  as  the 
second  item  on  the  first  record  (Appendix  A;  i = 1 , i=2)  of  the  parameter 
file  (*PARAMETERS. ) . 

3.  The  forest  unit  being  simulated  can  be  stratified  into  tree 
classes.  The  sample  data  can  be  organized  in  more  than  one  way;  a tree 
class  may  represent  one  tree,  a group  of  similar  trees,  or  a collection 
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of  similar  plots.  It  is  most  efficient  first  to  make  one  set  of  tree 
classes  that  reflect  average  conditions  over  the  forest  unit  being 
studied,  then  subidivide  by  plots  or  groups  of  plots  to  examine 
variability.  If  an  initial  condition  file  contains  more  tree  classes 
(records)  than  you  wish  to  use  in  a particular  simulation,  a subset  of 
the  file  can  be  selected  by  stipulating  the  first  and  last  tree-class 
number  (Figure  2,  item  9)  for  the  subsets  desired.  Up  to  10  subsets  can 
be  selected  (See  Appendix  A;  i =1 , j=3,...22). 

Now,  for  a specific  example,  assume  that  an  initial  condition  file 
contains  48  records,  numbered  consecutively.  The  first  24  records  are 
for  Douglas-fir  and  the  last  24  for  grand  fir.  You  want  to  simulate 
using  only  the  first  12  classes  from  each  host  species.  The  values  1, 
12,  25,  and  36  are  thus  entered  as  the  third  through  sixth  entries  on 
record  (card)  number  1.  All  remaining  entries  can  be  left  blank.  The 
second  entry  on  card  1 will  then  be  24,  which  is  the  number  of  tree 
classes  to  be  simulated.  If  all  48  tree  classes  are  desired,  then  the 
third  and  fourth  entries  on  the  first  record  should  be  1 and  48,  and  the 
second  entry  on  card  1 will  be  48. 

B.  Insect  Natural  Mortality  and  Fecundity 

As  the  model  insect  population  goes  through  its  annual  cycle, 
population  numbers  are  affected  by  mortality  from  natural  factors  and 
direct  control.  Control  will  be  discussed  in  Section  V-C.  Natural  lar- 
val mortality  is  divided  into  four  sources--background , predator/para- 
site, disease,  and  stress.  Background  mortality  is  the  same  for  each 
phase  and  is  the  assumed  phase  I net  rate.  Predator/parasite,  disease, 
and  stress  mortalities  are  assumed  to  act  independently,  both  of  each 
other  and  of  the  background  rate.  Predator/parasite  and  disease  rates 
are  entered  as  phase-  and  instar-specific  daily  rates.  Independence 
implies  the  following  multiplicative  rate,  mx,  for  a mortality  factor 
and  the  associated  survival,  sx: 

sx  = 1.0  - mx 

Net  daily  survival  (s)  is  the  product  of  the  individual  survival  factors 
(s  = sasb...  .).  Similarly,  survival  (S),  over  the  lOd  instar  is 
the  product  of  the  10  daily  survival  rates 

S = s s ^ • • ,sio* 

Table  5 gives  daily  and  instar  rates  from  each  factor,  as  well  as  the 
net  daily  rate  for  all  factors  other  than  stress  mortality.  Stress 
mortality  operates  only  when  the  population  is  forced  to  eat  old 
foliage.  This  occurs  after  all  new  foliage  has  been  consumed. 

Number  of  days  of  stress  mortality  must  be  known  (see  Section  V-G) 
to  incorporate  this  factor  in  net  instar  mortality.  For  example,  if 
five  days  of  feeding  on  foliage  occurs  in  the  fourth  instar  on  grand  fir 
during  phase  II: 
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Table  5.  (a)  daily  larval  mortality  rates  and  net  daily  rate  with  and 

without  stress  rate  included. 


Phase  or 

Instar 

Source 

host  species 

1 

2 

3 

4 

5 

6 

(a)  Daily 
background 

I,  II,  HI, 

IV 

0.02* 

0.02 

0.02 

0.02 

0.02 

0.02 

Predator  and 

I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

parasite 

II 

0.0 

0.0 

0.0 

0.0 

0.001 

0.001 

III 

0.001* 

0.002 

0.003 

0.010 

0.016 

0.042 

IV 

0.005 

0.006 

0.007 

0.021 

0.033 

0.056 

Di sease 

I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

II 

0.0 

0.0 

0.0 

0.0 

0.001 

0.001 

III 

0.002* 

0.003 

0.006 

0.013 

0.035 

0.028 

IV 

0.025 

0.028 

0.031 

0.034 

0.035 

0.035 

Food  stress 

Douglas-fi r 

0.92 

0.60 

0.07 

0.0 

0.0 

0.0 

grand  fir 

0.95 

0.70 

0.10 

0.02 

0.0 

0.0 

Net  daily 

I 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

rate  without 

II 

0.02 

0.02 

0.02 

0.02 

0.027 

0.027 

food  stress 

III 

0.023* 

0.025 

0.029 

0.042 

0.069 

0.087 

IV 

0.049 

0.053 

0.057 

0.073 

0.086 

0.107 

Net  daily 

I 

0.922 

0.608 

0.089 

0.020 

0.020 

0.020 

rate  including 

II 

0.922 

0.608 

0.089 

0.020 

0.022 

0.022 

Douglas-fir 

III 

0.922 

0.610 

0.097 

0.042 

0.069 

0.087 

stress 

IV 

0.924 

0.621 

0.123 

0.073 

0.086 

0.107 

Net  daily 

I 

0.951 

0.706 

0.118 

0.040 

0.020 

0.020 

rate  including 

II 

0.951 

0.706 

0.118 

0.040 

0.022 

0.022 

grand  fir 

III 

0.951 

0.707 

0.126 

0.062 

0.069 

0.087 

stress 

IV 

0.952 

0.716 

0.151 

0.092 

0.086 

0.107 

* (1.0-0.023)  = ( 1.0-0.02)  ( 1.0-0.001 ) ( 1.0-0.002 ) 
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Table  5.  (b)  instar  mortality  rates  and  the  net  instar  rates  excluding 

stress. 


Source 

Phase  or 
host  species 

Instar 

1 

2 

3 

4 

5 

6 

(b)  Instar 

Background 

I,  II,  III 

0.183 

0.183 

0.183 

0.183 

0.183 

0.183 

IV 

Predator  and 

I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

parasite 

II 

0.0 

0.0 

0.0 

0.0 

0.010 

0.010 

III 

0.010 

0.020 

0.030 

0.096 

0.149 

0.349 

IV 

0.049 

0.058 

0.068 

0.191 

0.285 

0.438 

Di sease 

I 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

II 

0.0 

0.0 

0.0 

0.0 

0.010 

0.010 

III 

0.020 

0.030 

0.058 

0.123 

0.300 

0.247 

IV 

0.224 

0.247 

0.270 

0.292 

0.300 

0.300 

Food  stress 

Douglas-fir 

l-io-ii 

1-10"4 

0.516 

0.0 

0.0 

0.0 

grand  fir 

1-10-13 

1-10-5 

0.651 

0.183 

0.0 

0.0 

Net  instar 

I 

0.183 

0.183 

0.183 

0.183 

0.183 

0.183 

rate  excluding 

II 

0.183 

0.183 

0.183 

0.183 

0.199 

0.199 

stress 

III 

0.207 

0.223 

0.253 

0.352 

0.513 

0.600 

mortal ity 

IV 

0.397 

0.421 

0.444 

0.532 

0.591 

0.678 
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S4net  ■ S4  W-02)\ 

= (0.817)  (0.98)b 

= (0.817)  (0.904) 

= (0.739) 

Here  S4  is  the  phase  II,  fourth-i nstar  survival,  excluding  stress. 

Thus  73.9  percent  of  the  number  starting  the  fourth  instar  will  finish 
it  (mortality  of  26.1  percent). 

Suppose  now  that  you  are  interested  in  the  rate  of  survival  (or 
mortality)  associated  with  disease  for  the  fourth  instar  of  phase  III. 
Because  daily  disease-related  mortality  rate  is  0.013  (Table  5a),  the 
(lOd)  instar  survival  rate  is  (1.0-0.013)10  = 0.98710  = o.877,  and 
the  (lOd)  instar  mortality  rate  is  1.0-0.877  = 0.123.  To  alter  the 
given  rate,  and  with  data  to  support  an  instar  mortality  rate  of  0.214 
(for  example,  21.4  percent  of  a rearing  sample  dies  of  disease  during 
the  fourth  instar),  the  daily  rate  would  be  obtained  for  entry  into 
*PARAMETERS.  as  follows: 


Sinstar=1*0-0-214 

Si nstar  = 0,786  ' S10d 

S.  = 0.786^  = 0.976  = 1.0-0.024, 
daily  * 

hence  the  daily  survival  rate  to  enter  in  replacement  of  0.013  in 
*PARAMETERS.  is  0.024. 


Occasion  mortalities  are  assumed  to  operate  similarly.  Only  one 
(natural)  rate  is  provided  for  each  occasion.  For  example,  all  natural 
mortality  factors  affecting  pupae  in  phase  II  are  assumed  to  total  62 
percent  (value  0.62  for  i =1 0 , j=2  in  Appendix  A).  If  pupal  mortality  is 
believed  to  be  different  during  phase  II  for  some  stand  (area)  being 
simulated,  the  second  entry  on  the  tenth  record  of  *PARAMETERS. 

(Appendix  A;  i = 1 0 , j=2)  can  be  replaced  with  the  new  value. 

Nominal  egg-mass  size  is  the  number  of  viable  eggs  produced  by  an 
adult  female  that  has  fed  for  the  full  larval  interval  on  new  foliage. 
This  number  is  reduced  by  one  percent  for  each  day  of  feeding  on  old 
foliage  and  further  reduced  if  the  model  branch  is  totally  defoliated. 
The  assumed  nominal  egg  mass  sizes  can  be  replaced  by  changing  the  last 
two  entries  of  record  10  (Appendix  A;  i = 1 0 , j=5,6)  or  the  first  two  on 
record  11  (Appendix  A;  i =1 1 , j = 1 , 2 ) . 

C.  Insect  Control  Mortality 

To  simulate  direct  control,  first  decide  the  phase  being  controlled, 
the  time  over  which  control  is  to  operate,  and  survival  for  the  affected 
stages.  Assume,  for  instance,  that  previously  accumulated  data  on  a 
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particular  pesticide  suggest  efficacy  of  97.4  percent.  The  population 
is  in  the  second  year  of  visible  (some  heavily)  defoliation.  Samples, 
taken  pre-  and  post-spray,  and  a subsample  indicate  that  about  half  the 
population  was  third  and  half  fourth  instars  before  spray,  and  half 
fourth  and  half  fifth  after  spray.  A daily  control  mortality  rate  is 
sought  that  will  give  97.4  percent  total  mortality  during  the  fourth 
instar  of  phase  III.  Because  fourth-i nstar  natural  mortality  parameters 
are  0.020,  0.013,  0.010  (Table  5): 


54th  instar 


10 

IT 

j=l 


Sd 


= <sd.> 
J 


10 


S4th  instar  = ( 1 - 0-0 .020 ) 1 °( 1 . 0-0 . 01 3 )10 ( 1 . 0-0. 01 0) 1 0 


(1.0-m  ) 

c4 


10 


And  from  the  efficacy  data:  S 
so  that: 


4th  instar 


= (1.0-0.974)  = 0.026, 


0.026  = (0. 980 )10( 0.987 )10( 0.990 )10(S  )10 

c4 


0.026  = (0.817)  (0.877)  (0.904)(s  ) 

c4 


10 


0.026  = (0.648) (s  ) 
c4 


10 


, vlO  0.026 

(s  ) = 

c4  0.648 


0.040 


l 

( s ) = (0.040)  10  = 0.775 

C4 


Thus  the  mortality  parameter  for  97.4  percent  control  in  the  fourth 
instar  of  phase  III  is  mc.  = 1.0-0.725  = 0.275,  or  27.5  percent  of 
the  daily  mortality  is  assumed  to  be  a direct  result  of  the  pesticide. 

If  control  acts  over  a longer  period,  you  must  decide  what  fraction 
of  the  control  mortality  should  be  included  in  each  instar  over  which 
the  control  is  assumed  to  act.  For  short-lived  chemical  agents,  you  may 
want  to  slow  control  in  earlier  instar(s)  by  assigning  the  highest  rate 
to  the  next  instar  after  the  simulated  application.  With  bacterial  or 
viral  control,  the  disease  will  spread--and  the  rate  might  increase  in 
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later  instars  (Stelzer  et  al . 1975,  Stelzer  et  al . 1977).  Daily  rates 
of  mortality  assigned  to  each  instar  should  reflect  existing  data  on  the 
efficacy  of  the  chemical  or  microbial  being  modeled. 


If,  in  the  above  example,  the  97.4  percent  efficacy  was  assumed  to 
operate  over  the  third  and  fourth  instars,  then: 

10  10 


(1.0-0.974)  = 0.026  = S3S4 


=( 1.0-0. 020) 10( 1.0-0. 006) 10( 1.0-0. 003) 10(l-m  ) 10 

C3 

• (1.0-0. 20) 10( 1.0-0. 013) 10( 1.0-0. 010) 10(l-m  ) 10 

c4 

= (0.817) (0.942) (0.970)  (1-m  )10 

c3 

• (0.817 ) (0.877) (0.904)  (1-m  )10 

c4 

0.026  = (0.4837 ) (1-m  )10(l-m  )10 

c3  c4 

i 

(1-m  ) (1-m  ) = (0.0537)  10  = 0.746. 

c3  c4 

When  defoliation  removes  all  new  foliage  before  the  end  of  the 
period  for  which  control  parameters  are  to  be  altered,  stress  mortality 
should  also  be  included  in  the  calculations. 

To  include  the  effects  of  control  on  pupal/adult  or  overwinter 
mortality,  the  same  multiplicative  rule  is  used: 

■*‘”mtotal " ^ *’mnatural  ^ ^ ^control  ^ * 

For  example,  if  the  data  indicate  pupal /adult  mortality  of  78  percent  as 
a result  of  phase  II  microbial  control,  and  phase  II  naturally  occurring 
pupal  mortality  is  0.620  (Table  3): 

(1.0-0.780)  = (1.0-0.620)  (l-mc), 

sol vi ng  for  mc  gi ves 

mc=  0.421. 

The  same  rule  is  used  for  combining  natural  overwinter  mortality  and 
control -i nduced  overwinter  mortality. 
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D.  Foliage  Nominal  Conditions 


Nominal  whole-crown  dry  weight  is  used  along  with  the  number  of 
trees  in  computing  stand  mean  insect  density  for  redistribution.  If  you 
want  mean  density  of  insects  per  unit  area,  then  crown  weights  and 
stocking  information  must  be  provided.  These  are  not  necessary  if 
redistribution  is  not  considered;  then,  the  second  and  third  entries  on 
each  tree  class  (Figure  2)  should  be  left  blank.  If  crown  weights  are 
needed,  then  whole-crown  foliage  weights  can  be  used;  they  are  provided 
by  the  following  regression  equations  (Brown  1978): 

d:  diameter  at  breast  height,  inches 

w:  foliage  dry  weight,  pounds 

Domi nant/codomi nant : 

grand  fir: 

w = exp  (1.3094  + 1.6076  ln(d) ) for  l<d<36 
1.592  + 0.059d 

w = 0.286  exp  (1.3094  + 1.6076  ln(d))  for  36<d<40. 

Douglas-f i r : 

w = 0.484  exp  (1.1368  + 1.5819  ln(d)-0.0210d)  for  l<d<17 

w = (0.4955d2-10.04)  exp(-0.0210d)  for  17<d<34 
Intermedi ate/suppressed : 

grand  fir  w = 1.5747  exp(1.6156  In  (d)-0.0544d) , l<d<12. 

Douglas-f i r w = 0.5977  exp(1.862  In (d)-0.0552d) , l<d<ll. 

The  Brown  publication  also  contains  foliage  information  for 
subalpine  fir.  Because  it  uses  only  data  from  the  northern  Rockies,  it 
does  not  contain  any  information  for  white  fir  or  species  in  other 
geographic  areas. 

Nominal  foliage  complements  for  the  model  branch  have  been  derived 
from  five  stands  each  from  Arizona,  New  Mexico,  and  the  northern  Rockies 
(Idaho  and  Montana),  as  well  as  10  stands  from  California  (Table  6). 
Complete  statistics  on  nominal  foliage  complements  for  Oregon  and 
Washington  have  not  been  fully  developed. 

Behavior  of  the  model  has  been  investigated  by  developing  an  array 
of  initial  conditions  that  covers  the  range  of  the  Douglas-fir  and  grand 
fir  data  and  the  range  of  phase  I initial  insect  densities  (Mason  1978; 
Hatch  and  Mika  1978).  From  this  array,  input-output  tables  have  been 
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Table  6.  The  means  and  standard  means  for  biomass  and  % new  foliage  for 
the  model  branch  by  geographic  regions  (data  entries  4 and  5 
for  the  initial  condition  file.  Fig.  2)  (from  Hatch  and  Mika 
1978). 


Host 

species 

National  Forest 
(State) 

Mean 

biomass 

S.E. 

(n)* 

Mean 
% new 

S.E. 

Douglas- 

Clearwater  (ID) 

219.2 

12.55 

(19) 

27.3 

1.39 

fir 

Flathead  (MT)** 

280.3 

15.83 

(20) 

40.4 

4.25 

Coeur  d'Alene  (ID) 

221.9 

16.54 

(10) 

26.6 

2.13 

Santa  Fe  (NM)*** 

221.9 

7.67 

(50) 

— 

— 

Tonto  (AZ) 

175.5 

7.18 

(50) 

21.4 

0.983 

Grand  fir 

Clearwater  (ID) 

235.7 

15.21 

(20) 

32.4 

1.58 

Coeur  d'Alene  (ID) 

209.8 

16.76 

(10) 

40.9 

1.08 

White  fir 

Eldorado  (CA) 

301.9 

15.77 

(30) 

19.9 

1.37 

Modoc  (CA) 

395.9 

12.39 

(50) 

19.8 

1.24 

Santa  Fe  (NM)*** 

303.1 

15.87 

(40) 

-- 

-- 

Stanislaus  (CA) 

208.2 

11.63 

(20) 

26.4 

2.00 

Tonto  (AZ) 

230.9 

7.21 

(50) 

23.2 

0.945 

* n:  The  number  of  samples.  A sample  consisted  of  three  18-inch 
branches,  one  outer  and  two  inner,  from  the  midbole  of  a tree  60 
feet  or  less  in  total  height. 

**  Previous  defoliation  of  portions  of  some  of  these  samples  may 
account  for  the  high  percent  new  foliage. 

***  Limited  current  year's  (1977)  foliage  was  observed  when  samples 

were  taken.  This  may  also  bias  the  total  foliage  biomass  for  this 
area. 
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produced  and  can  be  made  available  for  user  reference.  These  tables  use 
all  the  standard  default  values  given  in  Appendix  A. 

E.  Previous  Defoliation:  Estimating  Initial  Foliage  Conditions 

The  normal  procedure  for  using  the  model  has  been  to  simulate  the 
full  four-year  episode,  and  assume  an  initial  foliage  complement  equal 
to  the  nominal.  When  a simulation  is  to  start  in  a later  phase,  the 
user  must  locate  both  appropriate  nominal  foliage  conditions  and  initial 
insect  densities  for  earlier  years  that  will  match  collected  data^. 

From  this,  either  back  up  initial  conditions  to  phase  I or  use  the  ini- 
tial foliage  the  model  has  given  for  the  year  the  simulation  is  to 
begin. 

If  changes  are  made  in  the  file  that  would  affect  model  behavior 
before  the  point  where  an  accurate  simulation  is  desired,  the  above 
procedure  will  no  longer  apply.  The  user  should  then  make  an  initial 
set  of  simulations  to  assess  the  effects  of  changes  in  the  parameter 
file  and  choose  the  range  of  nominal  foliage  conditions  and  phase  I 
insect  densities  desired. 


F.  Insect  Population  Estimates  and  Model  Initial  Insect  Densities 

To  determine  the  initial  insect  density  that  will  give  a specific 
density  at  an  intermediate  point  within  any  one  year,  the  accumulated 
mortality  associated  with  the  given  period  must  be  used  and  the  derived 
equation  solved.  For  example,  if  the  data  give  an  average  of  52  third 
instars/1000  sq.  in.  in  phase  III  and  from  the  age  distribution  of  the 
sample,  this  number  is  assumed  to  correspond  to  the  midpoint  of  the 
third  instar,  then 


, 52  ^initial  SlS2s31s32s33s34s35’ 

where 


^initial : numbers  of  insects  (eggs), 

S^:  survival  for  first  instars, 

S^:  survival  for  second  instars, 

s3 1 "s 35 * c*a1^  survival  for  days  1-5  of  the  third 

i nstar 


From  Table  5,  S\  = 0.793,  S2  = 0.777,  and  S31  = ...  = S35  = 

0.9712.  Thus  the  initial  density  that  will  result  in  52  larvae  after 
five  days  of  the  third  instar  is  given  by: 


8 A set  of  behavior  tables  has  been  developed  and  is  available  through 
the  Methods  Application  Group. 
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52 


N 

i niti al 

(0.793)  (0.777)  (0.9712)u 

= 52  = 97.67 

(0.793)  (0.777)  (0.864) 

So  98  established  first  instars  (initial  eggs)  will  give  52  insects  half 
way  through  the  third  instar.  We  have  assumed  no  stress  mortality  here. 
If  stress  mortal ty  need  be  incorporated,  the  instar  rate  can  be  read 
from  Table  5,  or  stress  mortality  can  be  obtained  from  the  parameter 
file  (Appendix  A),  and  the  survivorship  equations  as  described  in  V-B 
and  V-C  above  can  be  used. 

G.  Details  of  a Simulation 


Table  4 may  be  produced  after  a simulation.  Using  two  data  files 
produced  during  a simulation  (*LUN10.  and  *LUN25.)  and  a specifications 
file  you  supply  (*SPECTAB4. ) , the  program  TABLE4  produces  desired 
details.  The  data  files  *LUN10.  and  *LUN25.  produced  during  a simu- 
lation are  normally  available  for  five  working  days  after  the  run  and 
may  be  preserved  (via  the  0SAVE  command)  for  a longer  period  if  desired 
(Section  III-C). 

Any  of  the  four  state  variables,  the  number  of  days  of  feeding  on 
old  foliage  or  the  number  of  days  spent  without  food,  are  available 
through  the  use  of  the  Table  4 program.  Figure  8 gives  the  input  codes 
that  must  be  specified  to  retrieve  this  information. 

To  obtain  any  of  these  details,  specify  the  total  number  of  tree 
classes  in  the  simulation,  the  number  of  items  desired,  and  the  ( i , j ) - 
code  for  each  item.  For  example,  if  24  tree  classes  had  been  used  in  a 
simulation  and  both  the  number  of  larvae  starting  the  third  instar  and 
new  foliage  biomass  at  the  end  of  feeding  are  of  interest,  then  input 
the  following: 

1.  24  (tree  classes) 

2.  2 (number  of  items  desired) 

3.  (4,3) 

4.  (8,1) 

The  format  for  these  specifications  and  the  necessary  file  assignments 
were  given  in  previous  sections. 
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Figure  8.  The  detailed  information  available  in  Table  4 and  the  specification  codes  to  obtain  it. 


i : occasion  code;  this 

j:  variable  code;  the 

number  is  one 

first  four  are 

greater  than  the 

the  state  variab' 

access  times 

es. 

of  Figure  2. 

Nw  j 

i N. 

New  foliage 

bi omass 

Old  foliage 

biomass 

Number  of 

insects 

Mean  larval 
biomass-^ 

No.  days  on 
old  foliage 

No.  days 

wi thout 

food 

Initiation 
of  year 

(1,1) 

(1,2) 

(1,3) 

2/ 

2/ 

1/ 

Start  first 
instar 

(2,1) 

(2,2) 

(2,3) 

(2,4) 

2/ 

2/ 

Start  second 
instar 

(3,1) 

(3,2) 

(3,3) 

(3,4) 

(3,5)—/ 

(3,6)—/ 

Start  third 
instar 

(4,1) 

(4,2) 

(4,3) 

(4,4) 

(4,5)—/ 

(4 ,6)—/ 

Start  fourth 
instar 

(5,1) 

(5,2) 

(5,3) 

(5,4) 

( 5,5)—/ 

( 5,6)—/ 

Start  fifth 
instar 

(6,1) 

(6,2) 

(6,3) 

(6,4) 

(6,5)—/ 

(6,6 )— 7 

Start  sixth 
instar 

(7,1) 

(7,2) 

(7,3) 

(7,4) 

( 7 , 5)—/ 

( 7 ,6)—/ 

End  larval 
stage,  start 
pupal 

(8,1) 

(8,2) 

(8,3) 

(8,4) 

(8,5)—/ 

(8,6)—/ 

Adult  females 
to  produce 
• eggs 

2/ 

2/ 

(9,3) 

2/ 

2/ 

2/ 

Before  over- 
winter 
mortal ity 

(10,1)^ 

( 10 ,2)— 7 

(10,3) 

2/ 

2/ 

2/ 

y Mean  larval  biomass  follows  a preassigned  trajectory  (according  to  initial  biomass  and  growth 
rates)  and  changes  only  when  the  model  branch  is  totally  defoliated,  i.e.,  item  (i,6). 


2/  The  areas  left  blank  indicate  no  usable  information  is  available. 
3/  New  foliage  (bud)  potential  for  the  following  year's  new  foliage. 
4/  Total  foliage  biomass  before  any  shedding  of  old  foliage. 

5/  Number  of  days  during  the  instar  just  completed. 
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APPENDIX  A 


The  parameters  of  *PARAMETERS.  and  their  default  values.  The 
indices  i and  j refer  to  the  record  or  card  number  (i)  and  the  item  (j) 
on  that  record. 


i 

j 

Val  ue 

Description 

Units 

1 

1 

1 

Beginning  phase  of  the  outbreak 
simulation 

yr 

1 

2 

** 

No.  of  tree  classes  to  be 
simulated 

none 

1 

3 

** 

First  tree  class  index  of  first 
subset 

none 

1 

4 

** 

Last  tree  class  index  of  first 
subset 

none 

1 

5 

— 

First  tree  class  index  of  second 
subset 

none 

1 6 

• 

• 

™ "" 

Last  tree  class  index  of  second 
subset 

none 

• 

1 

21 

-- 

First  tree  class  index  of  the 
tenth  subset 

none 

1 

22 

-- 

Last  tree  class  index  of  the 
tenth  subset 

**:  Option 

of  user 

but 

must  be  filled  in 

Option 

of  user 

and 

can  be  omitted 

Note: 

tree 

class  indices  must  be  less  than  1000 
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i 

2 

2 

2 

2 

2 

2 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 

4 

5 

5 

5 

5 

5 

5 

6 

6 

6 

6 

6 

6 

7 

7 

7 

7 

7 

7 


j Value  Description 


Units 


1 0.0 

2 0.0 

3 0.0 

4 0.0 

5 0.0 

6 0.0 


Instar 

1 

2 

3 

4 

5 

6 


Phase 


1 0.0  1 

2 0.0  2 

3 0.0  3 

4 0.0  4 

5 0.001  5 

6 0.001  6 


II 


Instar-specific 
daily  disease 
mortality  rate 


d-1 


1 0.002  1 

2 0.003  2 

3 0.006  3 

4 0.013  4 

5 0.035  5 

6 0.028  6 


1 0.025  1 

2 0.028  2 

3 0.031  3 

4 0.034  4 

5 0.035  5 

6 0.035  6 

1 0.0  1 

2 0.0  2 

3 0.0  3 

4 0.0  4 

5 0.0  5 

6 0.0  6 

1 0.0  1 

2 0.0  2 

3 0.0  3 

4 0.0  4 

5 0.001  5 

6 0.001  6 


Instar-speci f ic 
daily  predator/  d“l 
parasite  mortal  i ty 
rate 
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1 


Val  ue 


Description 


Units 


xx : 


8 1 0.001 

8 2 0.002 

8 3 0.003 

8 4 0.010 

8 5 0.016 

8 6 0.042 

9 1 0.005 

9 2 0.006 

9 3 0.007 

9 4 0.021 

9 5 0.033 

9 6 0.056 


Instar 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 


Phase 


III 


IV 


Instar-speci f ic 
daily  predator/  d~! 
parasite  mortal  ity 
rate 


10 

1 

0.50 

Phase 

I 

\ 

10 

2 

0.62 

II 

Natural  pupal 

10 

3 

0.75 

III 

> 

mortality  rate 

10 

4 

0.80 

IV 

10 

5 

200.0 

I 

'j 

Nomi nal 

10 

6 

200.0 

II 

/ 

Maximum  egg- 

11 

1 

150.0 

III 

mass  size; 

11 

2 

150.0 

IV 

) 

no.  of  eggs/mass 

occasion"! 


number 


11  3 xx 


blank  field 


none 


11 

4 

0.50 

I 

11 

5 

0.60 

II 

11 

6 

0.85 

III 

12 

1 

0.90 

IV 

Natural  occasion"! 

overwinter  mortality 

rate 


12 

12 

12 

12 

12 


0.0 

0.0 

0.0 

0.0 

0.0 


Redistribution  dispersal  factor 


I 

II 

III 

IV 


) 


Establ i shment 

mortality 

rate 


occasion 


-1 


location  is  not  used 
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j 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 


Val  ue 


Descri ption 


Units 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


Instar 

1 

2 

3 

4 

5 

6 


Phase 

\ 

I 


1 

2 

3 

4 

5 

6 


II 


Instar-speci f ic 
daily  control 
mortality  rates 


d-1 


1 

2 

3 

4 

5 

6 


III 


1 

2 

3 

4 

5 

6 


Phase 

0.0 

I ^ 

0.0 

11  f 

Control  pupal 

0.0 

III  J 

mortality  rate 

0.0 

IV  / 

occasion"! 


0.0  I \ 

0.0  II  1 Control  overwinter  occasion"! 

0.0  III  \ mortality  rate 

0.0  IV  ) 

xx 

XX 

XX 

XX 
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1 


J 


Val  lie 


Descri ption 


Units 


Instar 


19 

1 

0.92 

1 \ 

19 

2 

0.60 

2 \ 

Instar-specific  daily 

19 

3 

0.07 

3 

► stress  mortality  rate 

19 

4 

0.0 

4 

for  Douglas-fir 

19 

5 

0.0 

5 

19 

6 

0.0 

6 ) 

20 

1 

0.95 

20 

2 

0.70 

2 1 

Instar-specific  daily 

20 

3 

0.10 

3 

^ stress  mortality  rate 

20 

4 

0.02 

4 

for  grand  fir 

20 

5 

0.0 

5 

20 

6 

0.0 

6 i 

21 

1 

0.02 

1 'l 

! 

21 

2 

0.02 

2 

Instar-specific  daily 

21 

3 

0.02 

3 

^ background  mortality 

21 

4 

0.02 

4 

rate 

21 

5 

0.02 

5 1 

21 

6 

0.02 

6 J 

1 

22 

1 

5.40 

1 1 

1 

22 

2 

6.25 

2 

Instar-specific  destruction/ 

22 

3 

6.25 

3 

f consumption  ratio  for  new 

22 

4 

2.71 

4 

foliage 

22 

5 

2.27 

5 

1 

22 

6 

2.20 

6 J 

23 

1 

5.40 

1 ^ 

i 

23 

2 

6.25 

2 

Instar-specific  destruction/ 

23 

3 

6.25 

3 

l consumption  ratio  for  old 

23 

4 

3.69 

4 

[ foliage 

23 

5 

3.29 

5 

23 

6 

3.20 

6 

d-1 


d-1 


d-1 


none 


none 
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1 


J 


Val  ue 


Description 


Units 


24 

1 

1.19 

Assimilated  food  to  growth 
coeff ici ent 

none 

24 

2 

0.081 

Assimilated  food  to  respiration 
coefficient 

d-1 

24 

3 

0.10 

Assimilation  efficiency  factor 

none 

24 

4 

0.1147 

Instantaneous  growth  rate  for 
the  first  three  instars 

d-1 

24 

5 

0.0836 

Instantaneous  growth  rate  for 
the  fourth  instar 

d-1 

24 

6 

0.0625 

Instantaneous  growth  rate  for 
the  last  two  instars 

d-1 

25 

1 

0.197 

Initial  biomass  of  mean  larvae 
at  establishment 

mg 

25 

2 

52.063 

Final  biomass  of  mean  female 
larvae  at  pupation  if  all  food 
requirements  have  been  met 

mg 

25 

3 

0.02 

Coefficient  of  reduction  in  new 
foliage  potential  from  the  no. 
of  days  of  total  branch 
defol i at  ion 

d-1 

25 

4 

0.01 

Coefficient  of  reduction  in 
fecundity  from  no.  of  days  of 
feeding  on  old  foliage 

d-1 

25 

5 

XX 

25 

6 

0.50 

Proportion  of  larvae  that  are 
male  (and  pupate)  at  the  end  of 
the  fifth  instar 

none 
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APPENDIX  B 
Sample  Runstream 


A "+"  to  the  left  of  the  line  means  that  the  line  is  not  entered  by 
you  but  is  a response  from  the  computer.  The  prompt  ">"  from  the  com- 
puter indicates  that  it  is  ready  for  a line  of  input. 


UBQ123 

+ *UNIVAC  110  OPERATING  SYSTEM  VER.  33-R3-3F  (RSI)* 

>RUN  M01MAG, 1105205205  ,RUN1 
+ DATE:  052179  TIME:  145456 

+ THE  SUBJECT  FILE,  BROADCASTS,  HAS  BEEN  UPDATED  AS  OF  05/21/79  11:21:21 

+ READY 
>0QUAL  RUN1 
+ READY 

>@ADD,L  DFTM*RUNSTREAM. PRE1 
+ ODELETE  ,C  *LUN25. 

+ FURPUR  27R3A  E33  SL73R1  05/21/79  14:55:55 

+ *LUN25  IS  NOT  CATALOGUED  OR  ASSIGNED 

+ FAC  STATUS:  400010000000 

+ 0ASG,UP  *LUN25. 

+ READY 

+ OUSE  25,*LUN25 
+ READY 

+ ODELETE  ,C  *LUN10. 

+ FURPUR  27R3A  E33  SL73R1  05/21/79  14:56:10 

+ *LUN10  IS  NOT  CATALOGUED  OR  ASSIGNED 
+ FAC  STATUS:  400010000000 

+ OASG,UP  *LUN10. 

+ READY 

+ OUSE  10,*LUN10 
+ READY 

+ OASG,T  *1 
+ READY 


+ OASG ,T  *2 
+ READY 

+ OASG,T  *3 
+ READY 
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+ MSG  ,T  *4 
+ READY 

+ MSG,T  *5 
+ READY 

+ MSG  ,T  *6 
+ READY 

+ MSG,T  *7 
+ READY 

+ MSG,T  *8 
+ READY 

+ @ASG,T  *9 
+ READY 

+ @ASG ,T  *11 
+ READY 

+ MSG,T  *12 
+ READY 

+ MSG, UP  *15 
+ READY 

+ 0FREE  *15 
+ READY 

+ MSG ,A  *15 
+ READY 

+ @ASG,UP  *16 
+ READY 

+ @FREE  *16 
+ READY 

+ (MSG, A *16 
+ READY 

+ (MSG , UP  *17 
+ READY 

+ @FREE  *17 
+ READY 

+ MSG, A *17 
+ READY 
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+ @ASG ,T  *18 
+ READY 

+ @ASG,T  *19 
+ READY 

+ (PASGJ  *20 
+ READY 

+ @ASG, UP  *21 
+ READY 

+ (9FREE  *21 
+ READY 

+ @ASG,A  *21 
+ READY 

+ @ASG,T  *22 
+ READY 

+ @ASG,T  *23 
+ READY 

+ @ASG,T  *24 
+ READY 

+ @ASG ,A  DFTM*DFTMDF. 

+ FACILITY  WARNING  000000000200 

+ 0USE  13 ,DFTM*DFTMDF 
+ READY 

+ @ASG,A  DFTM*DFTMGF. 

+ READY 

+ @USE  14}DFTM*DFTMGF 
+ READY 

+ @ASG ,A  *IC. 

+ READY 

+ 0ASG ,A  *PARAMETERS. 

+ READY 

+ @USE  30 ,*IC 
+ READY 


+ 0USE  60,*PARAMETERS 
+ READY 

>0XQT  DFTM*0 UTB RE AK. NAMES 
>0ADD,L  DFTM*RUNSTREAM. P0ST2 

+ 0FREE  *2 
+ READY 


+ 0FREE  *3 
+ READY 


+ @FREE  *4 
+ READY 


+ 0FREE  *5 
+ READY 


+ 0FREE  *6 
+ READY 


+ 0FREE  *7 
+ READY 


+ 0FREE  *8 
+ READY 

+ 0FREE  *9 
+ READY 

+ 0FREE  *10 
+ READY 

+ 0FREE  *11 
+ READY 

+ 0FREE  *12 
+ READY 

+ 0FREE  *13 
+ READY 

+ 0FREE  *14 
+ READY 

+ 0FREE  *18 
+ READY 

+ 0FREE  *19 
+ READY 

+ 0FREE  *20 
+ READY 
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+ @FREE  *21 
+ READY 

+ @FREE  *22 
+ READY 

+ OFREE  *23 
+ READY 

+ @FREE  *24 
+ READY 


+ OFREE  *25 
+ READY 

>@SYM,U  *15. ,2,FCR104 
>@SYM,U  *16. ,2,FCR104 
>@SYM,U  *17. ,2,FCR104 
>@SYM,U  *IC.,2,FCR104 
>@SYM,U  *PARAMETERS. ,2,FCR104 
>@FIN 


RUNID:  M01MAG  ACCT:  1105205205 
TIME:  TOTAL:  00:01:48.219 

CAU:  00:00:05.563 

CC/ER:  00:01:03.068 


PROJECT:  RUN1 

I/O:  00:00:39.588 
WAIT:  00:04:43.977 


FCCC  RESOURCE  TIME  : 
IMAGES  READ:  76  PAGES: 
START:  14:54:56  MAY  21,  1979 


00:00:42.944 

26 

FIN:  15:06:19  MAY  21,1979 


SIZE:  01 2K 


EST  COST:  TOTAL:  00010.17 

CPU:  00002.22  I/O:  00007.20  CONNECT:  0000.75 

♦TERMINAL  INACTIVE* 

>@@TERM 


UBQ123 

+ *UNIVAC  1100  OPERATING  SYSTEM  VER.  33-R3-3F  (RSI)* 

>@RUN  M0MAG, 1105205205  ,RUN1 
+ DATE:  052179  TIME:  175304 

+ THE  SUBJECT  FILE,  BROADCASTS,  HAS  BEEN  UPDATED  AS  OF  05/21/79  14:55:50 

XPQUAL  RUN1 
+ READY 

>@ADD,L  DFTM*RUNSTREAM.TABLE4B 

+ @ASG ,A  *SPECTAB4. 

+ READY 


+ @ASG ,A  *LUN10. 
+ READY 
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+ (MSG, A *LUN25. 

+ READY 

+ ODELETE ,C  *TABLE4. 

+ FURPUR  27R3A  E33  SL73R1  05/21/79  17:54:27 

+ *TABLE4  IS  NOT  CATALOGUED  OR  ASSIGNED 

+ FAC  STATUS:  400010000000 

+ OASG,CP  *TABLE4. 

+ READY 

+ OUSE  11,*TABLE4 
+ READY 

+ OUSE  12,*SPECTAB4 
+ READY 

+ OUSE  10,*LUN10 
+ READY 

+ OUSE  25 ,*LUN25 
+ READY 

+ OXQT  DFTM*NTABLE4.A2 

+ THERE  WERE  15  TREE  CLASSES  COUNTED  IN  *LUN10 
+ REQUEST  IS  TO  PROCESS  THE  FIRST  15  TREE  CLASSES 
+ THERE  WERE  47  VALID  ITEMS  REQUESTED 

+ OFREE  11 
+ READY 

+ OFREE  *TABLE4. 

+ FACILITY  WARNING  1 00000000000 

+ OFREE  12 
+ READY 

+ OFREE  10 
+ READY 

+ 0FREE25 
+ READY 

>@SYM,U  *TABLE4. ,2,FCR104 
>OSYM,U  *SPECTAB4. ,2,FCR104 
>OFIN 
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PROJECT:  RUN1 


RUN  ID:  M01MAG  ACCT:  1105205205 
TIME:  TOTAL:  00:01:14.055 

CAU:  00:00:02.763  I/O:  00:00:51.947 

CC/ER:  00:00:19.344  WAIT:  00:03:17.562 


FCCC  RESOURCE  TIME  : 00:00:25.423 

IMAGES  READ:  21  PAGES:  12 

START:  17:53:04  MAY  21,1979  FIN:  17:58:45  MAY  21,1979 


SIZE:  01 1 K 

EST  COST:  TOTAL:  00010.84 


CPU:  00001.11  I/O:  00009.36  CONNECT:  0000.37 

^TERMINAL  INACTIVE* 

>@@TERM 
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APPENDIX  C 


Other  model  parameters.  Several  model  parameters  in  the  file 
*PARAMETERS.  already,  have  been  discussed.  Twenty-eight  other  model  para- 
meters have  not  been  addressed;  for  the  most  part,  these  are  parameters 
you  will  not  be  interested  in  manipulating.  You  should  know,  however, 
what  their  effects  are  on  the  model  and  simulation  output. 

Four  parameters  are  related  to  overwinter  mortality.  These  are  the 
establishment  mortality  rates.  At  present,  knowledge  of  changes  in 
insect  density  between  fall  oviposition  and  spring  establishment  of 
larvae  is  insufficient  to  quantify  separately  overwinter  reduction  and 
spring  establishment  loss.  The  total  is  incorporated  as  one  overwinter 
parameter.  This  overwinter  mortality  accounts  for  all  mortality  between 
oviposition  and  establishment  of  the  feeding  population  the  following 
summer. 

Redistribution  of  insects  between  tree  classes  is  obtained  by 
changing  the  redi stribution  dispersal  factor  from  its  default  value  of 
zero.  This  factor  indicates  what  proportion  of  the  variance  in  insect 
density  between  tree  classes  should  be  removed  each  year.  When  given 
the  value  0.25,  the  variation  between  tree  classes  (weighted  by  number 
of  trees  and  crown  weight  per  tree)  will  be  reduced  by  25  percent  for 
each  year  simulated,  A range  of  0.20  to  0.30  was  found  to  give 
reasonable  behavior  when  compared  with  a limited  data  set  (four  stands 
examined).  With  stand  inventory  and  insect  sample  data  sufficient  to 
utilize  this  model  feature,  0.25  is  the  suggested  value. 

As  larvae  develop,  a series  of  food  assimilation  parameters  are  used 
as  well  as  the  destruction/consumption  ratios  and  growth  rates  given  in 
Appendix  A.  The  derivation  of  the  assimilation  parameters  is  given  in 
Overton  et  al . (1978).  Destruction/consumption  ratios  are  derived  from 
Beckwi th— ("1978)  and  growth  rate  data  are  from  Mason  (personal 
communication).  The  growth  rates  are  instantaneous  daily  rates.  They 
are  used  to  determine  the  change  in  average  individual  larval  biomass 
(dry  weight).  The  population  is  assumed  to  consume  foliage  according  to 
the  number  of  insects,  the  size  of  the  average  insect,  and  the  daily 
growth  rate: 


kt(1) 

Consumption^  = Wt^.(e  -1)  A 
where 

: number  of  insects  feeding  on  day  t, 

Wt^  : dry  weight  of  the  average  insect  to  begin  t, 

kt(1) 

(e  -1)  : growth  increment  (dry  wt.)  proportion  to 

be  accumulated  on  one  day  t from  growth  rate  k^. 
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A 


: assimilation  factor  for  conversion  of  insect  dry  weight 
gain  into  foliage  dry  weight  consumption. 

Consumption. : total  dry  weight  consumed  by  the  insects 
in  maintaining  growth  op  day  t. 

The  assimilation  factor  A = 

J_  (k,  +^2)  ( 1 q ) 

k3  kt  103  mg 

Here: 

k t : assimilated  food  to  growth  coefficient  (Appendix  A;  i=24, 

1 J=D, 

k9:  assimilated  food  to  respiration  coefficient  (Appendix  A ; 

* 1=24,  j=2) 

k3:  assimilation  efficiency  (Appendix  A;  i=24,  j=3) 

k^:  instantaneous  growth  rate  for  day  t. 

Changes  in  the  assimilation  coefficients  or  growth  rates  should  not  be 
made  without  first  reviewing  the  procedure  by  which  these  were  derived 
(Overton  et  al . 1978).  Also,  changes  in  growth  rates  should  not  be  made 
without  also  changing  the  initial  or  final  mean  larval  biomass 
accordi ngly. 

Consumption  is  transformed  to  foliage  destroyed  daily  by  use  of  the 
destruction/consumption  ratio  (Appendix  A;  i = 2 2 , 23).  This  ratio  is 
instar-  and  foliage  age-  (new,  old)  specific  but  not  host  species 
speci fic. 

Half  the  insects  are  assumed  to  be  males  and  pupate  at  the  end  of 
the  fifth  instar  (Appendix  A;  i =25 , j=6).  Assuming  that  male  and  female 
pupae  are  subject  to  the  same  mortality  factors,  females  will  be  sub- 
jected to  further  mortality  in  the  additional  i nstar--hence  differences 
in  adult  sex  ratio. 

The  remaining  two  parameters  are  coefficients  for  reducing  the 
ability  of  the  insect  and  host  to  reproduce.  First,  new  foliage  poten- 
tial, the  ability  of  a tree  to  flush  buds  set  in  the  previous  year,  is 
reduced  by  two  percent  (Appendix  A;  i=25,  j=3)  for  each  day  the  larvae 
would  be  actively  feeding  past  the  time  of  total  defoliation  (of  the 
model  branch).  This  accounts  for  buds  destroyed  by  larvae  during  this 
period  as  well  as  the  buds  killed  indirectly  by  stress  from  early  total 
defoliation.  New  foliage  potential  is  also  reduced  as  a result  of  defo- 
liation when  the  model  branch  is  only  partially  defoliated. 
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Second,  nominal  fecundity  is  reduced  by  one  percent  (Appendix  A; 
i=25,  j=4)  for  each  day  the  population  feeds  on  old  foliage.  This 
accounts  for  the  indirect  effects  of  food  quality  on  larval  development 
and  consequent  reduced  egg-mass  size. 
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APPENDIX  D 


The  two  host  effects  files,  DFTM*DFTMDF.  and  DFTM*DFTMGF. , are  read 
in  from  logical  units  13  and  14,  respecti vely,  using  the  following  two 
formats:  the  first  seven  records  use  (8F6.3);  the  last  two  records  use 
(6F6.3) . Figure  9 shows  the  standard  DFTM*DFTMDF.  and  DFTM*DFTMGF. 
files.  For  each  host  species,  the  effects  of  defoliation  were  deter- 
mined by  analysis  of  impact  plot  data  to  produce  an  array,  A,  of  expec- 
tations of  the  various  measurements  taken  as  follows: 


A ( i , j ) i-row:  card  or  record  number  in  the  file, 

j-column:  entry  number  across  a record. 

For  i=l,...,7  the  following  correspondences  exist: 

1:  all  parameters  correspond  to  trees  (a  tree  class)  with  less  than 
15  percent  of  the  crown  totally  defoliated, 

2:  all  parameters  correspond  to  trees  with  between  15  and  35 
percent  of  the  crown  totally  defoliated, 

3:  parameters  for  trees  with  between  35  and  65  percent  of  the  crown 
totally  defoliated, 

4:  parameters  for  trees  with  between  65  and  85  percent  of  the  crown 
total ly  defol i ated , 

5:  parameters  for  trees  with  between  85  and  95  percent  of  the  crown 
totally  defoliated, 

6:  parameters  for  trees  with  between  95  and  99.5  percent  of  the 
crown  totally  defoliated, 

7:  parameters  for  trees  with  between  99.5  and  100  percent  of  the 
crown  totally  defoliated, 

and  for  each  of  the  above  values  of  i,  the  following  definitions  of 

A(i,j)  for  j=l,...,8  exists: 

1:  expected  proportion  of  trees  for  which  mortality  is  directly 
attributable  to  defoliation, 

2:  expected  proportion  of  trees  with  no  top-kill  damage, 

3:  expected  proportion  of  trees  with  leader-kill,  one  year's  growth 

only 

4:  expected  proportion  of  trees  with  top-kill  less  than  10  percent, 

5:  expected  proportion  of  trees  with  top-kill  less  than  25  percent, 
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Figure  9.  (a)  the  standard  DFTM*DFTMDF.  file 
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0 
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0.  460 

0. 
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0. 

08 1 
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0. 
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0 

351 

0 

6 1 9 

0 

923 

0.  460 

0 

1 89 

0 

081 
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o. 
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0 

351 

0 
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o 

014 

0.  005 

o. 
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0. 

030 

0.  053 

0 
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o 

029 

0 021 

0. 

032 

0. 

000 

0.  000 

o. 

143 

(b) 

the 

standard 

DFTM*DFTMGF. 

■file. 

0 

000 

0. 

879 

0. 

073 

0 048 

0.  048 

0. 

048 

0. 

048 

0. 

892 

o 

000 

0 

762 

0. 

143 

0.  072 

0.  072 

0. 

072 

0. 

095 

0. 

80 1 

o 

009 

0 

555 

0 

191 

0 222 

0.  222 

0 

254 

0. 

254 

0 

754 

o 

0 2 c' 

0 

457 

0. 

257 

0.  171 

0.  171 

0. 

286 

0. 

286 

0. 

727 

o 

173 

0. 

475 

o. 

175 

0 075 

0.  075 

0 

200 

0 

350 

0 

698 

o 

477 

o 

462 

o. 

308 

0.  076 

0 076 

0 

076 

0. 

230 

0. 

572 

o. 

925 

0. 

462 

o. 

308 

0.  076 

0.  076 

0. 

076 

0. 

230 

0. 

572 

o 

068 

o. 

000 

0. 

034 

0.  333 

0 250 

0 

667 

o 

076 

0 

049 

o. 

000 

0.  145 

0.  000 

0. 

667 
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6:  expected  proportion  of  trees  with  top-kill  less  than  50  percent, 
7:  expected  proportion  of  trees  with  some  top-kill, 

8:  expected  loss  in  diameter  growth  from  defoliation. 


The  last  two  records  in  these  files  contain  the  information 
necessary  to  remove  secondary  mortality  from  the  arrays  described  above. 
They  are  as  follows: 


A(i,j):  i =8 : 

expectations  of  secondary  mortality  associated  with 
bark  beetle  secondary  attack. 

9: 

expectations  of  secondary  mortality  from  all  other 
causes,  e.g.,  windthrow,  winterkill,  natural 
attrition. 

j=l: 

expectation  for  trees  without  top-kill, 

2: 

expectation  for  trees  with  leader-kill  only. 

3: 

expectation  for  trees  with  more  than  leader-kill,  but 
less  than  10  percent  top  kill. 

4: 

expectation  for  trees  with  top-kill  between  10  and 

25  percent, 

5: 

expectation  for  trees  with  top-kill  between  25  and 

50  percent. 

6: 

expectation  for  trees  with  top-kill  greater  than  50 
percent. 

Mortality  and  degree  of  top-kill  and  growth  loss  have  been 
incorporated  as  the  effects  of  varying  ranges  of  defoliation.  The  model 
follows  the  population  of  insects  through  four  years  on  the  model  branch 
and  then  translates  model  branch  defoliation  into  percent  crown  totally 
defoliated.  From  this  variable,  one  of  the  host-specific  vectors  of 

effects  is  obtained  from  the  inputs  given  in  DFTM*DFTMDF.  or 
DFTM*DFTMGF.  (Section  II-B). 

Direct  mortality  has  been  measured  by  tallying  all  trees  dead  the 
year  after  peak  defoliation  that  show  no  signs  of  secondary  beetle 
attacks.  Any  tree  showing  signs  of  beetle  galleries  was  assumed  to  be 
killed  by  bark  beetles.  Trees  were  followed  for  an  additional  four 
years  after  tussock  moth  population  collapse.  All  trees  subsequently 
dying  and  containing  beetle  galleries  at  death  were  assumed  killed  by 
beetles.  All  other  mortalities  in  these  years  were  incorporated  into 
the  "other  causes"  category  (Wickman  1977). 
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The  two  forms  of  secondary  mortality  are  assumed  to  act 
independently,  hence  the  two  vectors  are  combined  by  the  rule: 

(mi,j  + m2,j  fY 

The  scalar-  or  dot-product  of  this  vector  with  the  input  vector  of 
residuals  of  direct  mortality  is  then  used  to  compute  secondary 
mortality  for  each  defoliation  class: 


Yi  J+l 


,7. 


Here  r- j j is  the  vector  of  residuals  (ith  record  of  DFTM*DFTMGF.  or 
DFTM*DFmDF. ) and  Sj  is  the  derived  secondary  mortality  that  will  be 
used  in  Table  2 of  simulation  output. 


After  secondary  mortality  has  been  subtracted  from  the  residuals  of 
direct  mortality,  the  remaining  numbers  are  the  expected  proportions 
that  are  to  receive  various  levels  of  top-kill,  and  one  class  that  indi- 
cates the  proportion  receiving  only  growth  reductions. 

Growth  reductions  attributable  to  defoliation  are  changes  in 
diameter  growth  at  breast  height  and  height  growth.  At  present,  only 
diameter  growth  has  been  sampled.  Until  further  data  become  available, 
the  reduction  multiplier  for  normal  diameter  growth  will  be  used  for 
reducing  height  growth.  These  multipliers  are  measures  of  the  propor- 
tion of  nominal  growth  attained  after  the  given  range  of  defoliation. 


The  parameters  in  either  of  these  files  can  be  altered 
to  affect  any  measured  or  theorized  changes  in  defoliation 


by  the  users 
effects. 
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